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ABSTRACT 

We report on spectroscopic analyses of three K giants previously suggested 
to be Li-rich: IRAS 12327-6523, IRAS 13539-4153, and IRAS 17596-3952. High- 
resolution optical spectra and the LTE model atmospheres are used to derive 
the stellar parameters: {Tes, log g, [Fe/H]), elemental abundances, and the iso¬ 
topic ratio IRAS 13539-4153 shows an extremely high Li abundance of 

loge(Li) 4.2, a value ten times more than the present Li abundance in the local 
interstellar medium. This is the third highest Li abundance yet reported for a K 
giant. IRAS 12327-6523 shows a Li abundances of Ioge(Li)f» 1.4. IRAS 17596- 
3952 is a rapidly rotating (Rsinz 35 km s“^) K giant with loge(Li) 2.2. 
Infrared photometry which shows the presence of an IR excess suggesting mass- 
loss. A comparison is made between these three stars and previously recognized 
Li-rich giants. 

Subject headings: stars:abundances-stars: Li-rich-stars: carbon isotopic ratios 
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1. Introduction 

The maximum lithium abuudauce of maiu sequeuce stars of uear-solar metallicity spaus 
a small rauge: loge(Li) 3.3 for [Fe/H] ^ 0.0 to 3.1 for [Fe/H] ^ —0.3 (Lambert & Reddy 
2003). This lithium abuudauce is ideutihed with that of the iuterstellar gas from which 
the stars formed. Lithium is predicted to be destroyed throughout a maiu sequeuce star 
except iu the outermost layers (1% - 2% by mass). Eveu iu this thiu skiu which iucludes the 
atmosphere, lithium may be depleted by processes uot yet fully uuderstood. The Suu, for 
example has a lithium abuudauce about a factor of 100 less thau the meteoritic abuudauce. 

Wheu the maiu sequeuce star evolves to become a red giaut, the couvective euvelope 
grows aud dilutes the lithium such that its surface abuudauce is greatly reduced. Ibeu 
(1967a,b) predicted that giauts of approximately solar metallicity would experieuce dilutiou 
by a factor of 1.8 dex at 3Mq aud 1.5 dex at IMq. This dilutiou occurs at the ouset of the first 
dredge-up resultiug from growth of the couvective euvelope which also results iu predicted 
chauges of the surface (elemeutal aud isotopic) abuudauces of helium, beryllium, borou, 
carbou, uitrogeu, aud oxygeu. Adoptiug the predicted lithium dilutiou aud the maximum 
observed lithium abuudauce of maiu sequeuce stars, lithium iu red giauts is expected uot to 
exceed loge(Li) Ri 1.5. A red giaut star that has evolved through the first dredge-up but has 
a surface abuudauce greater thau loge(Li) r^ 1.5 is declared to be ‘Li-rich’. 

Lithium-rich K giauts are rare. Siuce the discovery of the hrst Li-rich K giaut, HD 112127, 
two decades ago (Wallersteiu & Suedeu 1982), the collectiou of Li-rich stars has growu to 
ouly about three dozeu. Mauy like HD 112127 were discovered sereudipitously. Browu et 
ah (1989) uudertook a systematic search for Li-rich stars amoug K aud late G-type giauts. 
They fouud at most teu Li-rich giauts out of the 670 stars observed. Auother importaut 
survey kuowu as “Pico Dias Survey” (PDS) by Gregorio-Hetem et al. (1992) fouud, as a 
byproduct, poteutially Li-rich giauts (de la Reza et al. 1997). A few of these PDS caudidate 
Li-rich K giauts have beeu spectroscopically studied to determiue their lithium abuudauce 
(e.g; Reddy et al. 2002a; Drake et al. 2002). 

Here, we report results from the hrst high-resolutiou spectroscopic aualyses of three 
Li-rich caudidates takeu from de la Reza et al. (1997) who proposed the stars as probably 
Li-rich from the observatiou that the streugth of the Lil 6707 A doublet was equal to or 
greater thau the streugth of the Gal 6718 A hue. It turus out that the lithium abuudauces 
of the three stars are quite differeut. The most extreme example is IRAS 13539-4153 (also 
kuowu as PDS 68) with loge(Li) Rid.l, au abuudauce almost 1 dex higher thau the star’s 
likely iuitial lithium abuudauce. The secoud star IRAS 17596-3952 has loge(Li) r^ 2.2 aud is 
most probably a Li-rich giaut, i.e., a hrst dredge-up giaut with lithium iu excess of the limit 
loge(Li) ^ 1.5. The third star IRAS 12327-6523 has loge(Li) Ri 1.4 aud so is uot a certaiu 
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Li-rich star; it is possibly a post-first dredge-up star with lithium depleted to the predicted 
level. Following presentation and discussion of the abundance analysis, we comment on the 
similarities between the new and previously known Li-rich giants. 


2. Observations 

Spectra of the stars listed in Table 1 were obtained with the 4 m Blanco telescope and 
the Cassegrain echelle spectograph at the Cerro Tololo Inter-American Observatory in Chile. 
In addition to obtaining spectra of the three Li-rich candidates, we obtained a high-quality 
spectrum of HD 19745, a Li-rich giant analysed previously by de la Reza & da Silva (1995). 

The spectrograph was set to cover the wavelength interval 5000 A to 8200 A without 
gaps at a resolving power of i? ~ 35,000. Each star was exposed thrice with a total integration 
time ranging from 60 min to 90 min. Raw two-dimensional spectra were reduced to one¬ 
dimensional wavelength-calibrated spectra using appropriate reduction tools available in the 
I RAF package. A Th-Ar hollow cathode spectrum provided the wavelength calibration. In 
each spectrum of a Li-rich candidate star, the signal-to-noise at 6707 A is between 100 and 
200 . 


3. Analysis 
3.1. Stellar Properties 

Stellar parameters Teg, logg, and [Fe/H]) were obtained from the spectra by a standard 
procedure using ATLAS model atmospheres computed using the no convective overshoot 
option (Kurucz 1994). The effective temperature (Tefr) was obtained from the requirement 
that all Fel lines give the same abundance independent of the lower state’s excitation po¬ 
tential (LEP). About 40 Fel lines spanning the range LEPs of 1 to 5 eV were measured for 
IRAS 12327-6523 and IRAS 13539-4153 but fewer lines were measureable for the rapidly- 
rotating K giant IRAS 17596-3952. Lines are weak to moderately strong {W\ < 150 mA). 
Surface gravity (logg) was obtained with the condition that the Fe l and Fe ll lines give the 
same abundance. The microturbulent velocity (.^t) was determined by requiring that the 
Fel lines and Nil lines give the same (Fe or Ni) abundance irrespective of their equivalent 
width. The parameters were determined by an iterative process using models selected from 
the ATLAS grid. Adopted parameters are given in Table 1. Our results for HD 19745 are 
in good agreement with the previous determination by Castilho et al. (2000). 
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A full abundance analysis was made with the selected model atmospheres. Lithium 
abundance and the ratio are discussed below. Seventeen elements from Na to Nd 

in addition to Fe were included in the analysis. Given that [Fe/H] ~ 0 for the stars which 
probably belong to the thin disk, one expects [X/Fe] ~ 0 for the 17 elements. This is what 
is fonnd to within the measurement errors, except for the suggestion that [Ca/Fe] ~ —0.3 
for IRAS 12327-6523 and IRAS 13539-4153. The detailed analysis was not attempted for 
the rapidly rotating IRAS 17596-3952. 

Inspection of the spectra shows two nnusual aspects worthy of comment: unusnally 
broad lines of IRAS 17596-3952, and an asymmetric Ha prohle of IRAS 12327-6523. 

The lines of IRAS 17596-3952 are obviously broader than those of the other two stars 
which have line widths typical of K giants. The additional broadening is attributed to a 
more rapid than average rotational velocity. To determine the projected rotational veloc¬ 
ity (Rsini), we match the observed spectrum to a synthetic spectrnm computed with the 
following broadening parameters: the velocity Rsini, a macrotnrbulent velocity Vt, and the 
instrumental prohle width. The value Vt = 3 km s“^ was adopted, as found to be typical 
of K giants (Gray 1989). The instrnmental prohle width was taken as Gaussian with a 
FWHM = 0.21 A at 6707 A as measnred from the Th emission lines in the hollow cathode 
spectrnm. The non-Gaussian form (Gray 1976) of the broadening resulting from solid body 
stellar rotation is incorporated in the line formation code MOOG (Sneden 1973). Unblended 
Fel lines at 6703.5 A, 6705.1 A and 6170 A were selected for the estimation of Usini. A 
X^-test was used to select the best value of Usini. The mean value from the hts to the three 
Fe I lines is given in Table 1 for each star. 

Prohles of the Ha line in the four stars and Arctnrus are shown in Figure 1. It is clear 
that IRAS 12327-6523’s Ha line is distinctly asymmetric with the line core shifted to the 
bine relative to the photospheric lines. An asymmetric Ha line is seen in some other Li-rich 
giants (Drake et al. 2002). 


3.2. Lithium Abundance 

An immediate indication that the stars might be Li-rich was fonnd at the telescope from 
the obvious presence of a strong absorption line attribntable to the Li l resonance doublet at 
6707.8 A. Gloser inspection reveals also the presence of the excited Li l line at 6103.6 A in the 
spectra of IRAS 13539-4153 and HD 19745. Lithinm abnndances were fonnd by spectrum 
synthesis. 

Synthesis of the 6707 A feature used the critically examined line list in the vicinity of 
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the resonance line taken from Reddy et al. (2002b) with the component structure of the Li l 
resonance doublet taken from Hobbs, Thorburn, & Rebull (1999). Lithium is assumed to be 
purely ^Li. The line list for 10 A centered on the Li l line A 6707 A was tested by successfully 
reproducing the ultra-high [R ^ 150,000) resolution spectrum (Hinkle et al. 2000) of the 
very Li-poor K-giant Arcturus (see for details: Reddy et al. (2002a)). Then, synthetic 
spectra were computed and matched to the spectrum of each program star by adjusting 
the lithium abundance to £t the 6707.8 A feature. Comparisons of synthetic and observed 
spectra are shown in Figure 2 & 3. Lithium abundances derived from the LTE analysis are 
given in Table 2. Uncertainties in the derived lithium abundances (Table 2) are estimated 
from uncertainties in the atmospheric parameters: ST^s = 100 K, (51ogg = 0.5, S^t = 0.5 km 
s“\ and (5[M/H] = 0.25 dex lead uncertainties of 0.25, 0.05, 0.05, and 0.03 dex, respectively. 

For IRAS 13539-4153 and HD 19745, the lithium abundance has been obtained also 
using the subordinate Lil line at 6103 A. The line list for this region was compiled using 
Kurucz’s line list (Kurucz 1994). In Figure 4, we compare synthetic and observed spectra. 
The lithium abundances (Table 2) derived from the 6707 A and 6103 A lines are in good 
agreement for HD 19745, and in good agreement too with previously published results by 
de la Reza & da Silva (1995). 

In Figure 2, we draw attention to a failure to fit the Lil 6707 A profile of IRAS 13539- 
4153. The line is broader for its depth than calculated: the abundance loge(Li) = 2.35 
fits the core but loge(Li) = 4.1 is required to fit the wings, but this abundance produces 
a central depth greater than observed. The 6103A Lil line is matched by the abundance 
loge(Li) = 4.1 and the Usinf of the metal lines. We suppose that this latter abundance is the 
stellar abundance and that the 6707 A line’s core is affected by emission due to chromospheric 
activity or by non-LTE effects. It is unclear why IRAS 13539-4153’s 6707 A line is so affected 
but HD 19745’s line of comparable strength is not. (IRAS 13539-4153’s 6707 A profile is 
well fit by increasing the Usini to 16 km s“^ but this broadening fits no other lines in the 
spectrum (see Figure 2).) 

An additional uncertainty arises from our assumption of LTE for line formation. For 
the super Li-rich giants HD 19745 and IRAS 13539-4151 non-LTE corrections indicate a 
correction of 0.1 dex upward for the 6103 A line and 0.2 dex downward for the 6707 A line 
(Carlsson et al. 1994). For the less lithium-rich stars IRAS 12327-6523 and IRAS 17596- 
3952, the non-LTE correction is around 0.1 dex upward for the 6707 A line. Since these 
corrections are small and subject to their uncertainties, we include them as part of the 
overall abundance uncertainty. (The non-LTE abundances are also given in Table 2.) The 
total error given in Table 2 is the quadratic sum of all the errors discussed above. 
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3.3. The 12C/13C ratio 

Standard models make predictions about the first dredge-up’s effects on the surface 
abundances of Li, C, N, and O. Observations of normal red giants conhrm these predictions, 
qualitatively, if not always quantitatively. Interpretation of lithium abundances is assisted 
and possibly complicated by knowledge of the C, N, and O abundances. Here, we restrict 
the analysis to the carbon isotopic ratio One expects an unevolved star to have a 

ratio close to ~ 90, the solar photospheric value. A lower ratio (~ 20) is predicted 

for giants after the hrst dredge-up. Some giants show even lower ratio. 

The ratio was extracted from a complex of lines at 8005.7 A and adjacent 

lines (Figure 5). The ratio is determined by spectrum synthesis assuming 

solar abundances (Fodders 2003) for C and N. For HD 19745 we obtained = 16±2 

which is in good agreement with the ratio of 15 derived by da Silva, de la Reza, & Barbuy 
(1995). For IRAS 12327-6523, we obtain = 6±1. A limit of > 20 was found 

for IRAS 13539-4153. Spectral lines of IRAS 17596-3952 are too smeared, due to its high 
rotation, to estimate the ratio. The low ratios conhrm that the Li-rich giants have 

experienced the hrst dredge-up. 


4. The Li-rich clans 

When considered with the population of held K giants, our trio stand apart by reason 
of their high lithium abundance. That is why they were selected for an abundance analysis. 
The question we consider here is how the trio ht in with the previously analyzed Li-rich 
giants. 

Charbonnel & Balachandran (2000) took advantage of Hipparcos parallaxes to estimate 
the absolute luminosity of Li-rich and Li-normal red giants. They suggested that the Li-rich 
giants could be divided into three clans: warm giants for which lithium dilution has not 
been completed; giants at the bump in the luminosity function; giants in the early phases of 
asymptotic giant branch (AGB) evolution - see Figure 6. 

Our trio (and HD 19745) are not in the Hipparcos catalogue. Therefore, we estimate their 
luminosities using the spectroscopically-derived logg and Tes and computed evolutionary 
tracks for different stellar masses and metallicities (Girardi et al. 2000). By interpolation, we 
found the track (i.e., stellar mass) for which the logg and Tgfj best agreed with the observed 
log g and Tgfj. The luminosity is then read from the track. Estimated luminosities and 
masses are given in Table 1. The estimated errors due to uncertainties in model parameters 
(hTgfj = 100 K, hlogg = 0.25, and 6[Fe/H] = 0.2) are shown in Table 1. According to 
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their luminosity and effective temperature the three stars and HD 19745 have completed the 
dilution associated with the hrst dredge-up. The stars are shown in Figure 6 with others 
from our earlier studies. For stars taken from earlier studies uncertainties of ST^s = 100 K 
and (51ogg = 0.25 are assumed (Figure 6). 

The two most Li-rich giants IRAS 13539-4153 and IRAS 17596-3952 belong to the 
‘bump’ clan. By the spectroscopically estimated luminosity, we place HD 19745 also among 
this clan. The range of lithium abundances, ratios, and Rsini covered by these three 

stars are spanned by Charbonnel & Balachandran’s six stars: IRAS 17596-3952’s Rsini with 
35 km s“^ extends the upper limit previously set at 25 km s“^. 

IRAS12327-6523 is assigned to Charbonnel & Balachandran’s ‘early-AGB’ clan with 
its original hve members with lithium abundances ranging from loge(Li) = 1.5 to 2.8. 
IRAS 12327-6523 is now the least lithium rich of the clan. Strictly speaking, the star is 
not lithium rich because standard dilution by the first dredge-up could have left this rela¬ 
tively massive star with the now observed abundance. However, standard models of the first 
dredge-up cannot account for the low (=6) ratio; a ratio of 20 - 30 is predicted. It 

remains to be shown that the lower than expected carbon isotopic ratio can be achieved by 
modifications to standard models without severe reduction of the surface lithium abundance. 

Gregorio-Hetem et ah (1992, also Gregorio-Hetem, Gastilho, & Barbuy 1993) found 
that the great majority of the Li-rich giants possess an infrared excess. Drake et al. (2002) 
noted that every Li-rich rapid rotator in a sample of 20 had an infrared excess. Far-IR fluxes 
from the IRAS catalog show that our three stars (and HD 19745) have a strong excess at 12 
yum and longer wavelengths, as expected from their designation as an IRAS source. Fluxes 
estimated from the 2MASS J, H, and K magnitudes show that the IR excess does not extend 
to these shorter wavelengths. Our three Li-rich stars, therefore, continue the correlation 
between Li-richness and IR excess, and, in particular, IRAS 17596-3952 strengthens the 
impression that all Li-rich rapid rotators have a far IR excess. 


5. Concluding remarks 

We have performed the first detailed analysis of three Li-rich candidate K giants and 
a reanalysis of HD 19745. Two candidates are confirmed as certainly Li-rich: IRAS 13539- 
4153 is among the most Li-rich giants yet analysed. The status of IRAS 12327-6523 as 
Li-rich might be questioned. Three of the four giants - HD 19745, IRAS 13539-4153, and 
IRAS 17596-3952 - belong to RGB-clump collection of Li-rich giants and IRAS 12327-6523 
belongs to the early-AGB collection in the H-R diagram. This new trio and the other Li-rich 



giants analyzed previously by us confirm that characteristics of the RGB-bump clan that set 
them apart from other red giants: lithium abundance as high as log{Li) ~ 4, a far-infrared 
excess, and in the mean a high rotational velocity. 

Charbonnel & Balachandran (2000) linked the RGB-bump clan to the RGB bump where 
evolution along the giant branch is slowed when the H-burning shell burns through the 
molecular weight discontinuity left from the initial growth of the convective envelope which 
provided the first dredge-up. (The Li-rich giants at higher luminosity and of higher mass were 
associated with a molecular weight discontinuity predicted to occur in early-AGB stars as a 
result of the second dredge-up.). Erasure of the molecular weight gradient makes it possible 
for ‘extra’ mixing to occur between the convective envelope and the top of the H-burning 
shell. Mixing results in the conversion of ^He to ^Be and then ^Li by the Cameron-Fowler 
(1971) mechanism. This extra mixing may be internally or externally triggered. A proposed 
trigger should explain the concentration in luminosity of the Li-rich giants, the tendency for 
Li-rich giants to be rapid-rotators and the association with a far infrared-excess. A reader 
challenged to take up the search for the trigger might profitably read the proposals advanced 
by Palacios, Charbonnel, & Forestini (2001) and Denissenkov & Herwig (2004). 


David L. Lambert acknowledges the support of the Robert A. Welch Foundation of 
Houston, Texas. We would like to acknowledge with thanks use of 2MASS and NASA ADS 
data bases. 
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Table 1. 


Properties of the four Li-rich K giants. 


Star 

'Teft 

K 

log 9 

cm 

km 

hkelio 

km 

[Fe/H]“ 

Vsini 

km 

Wa(A6707) 

mA) 

M/M© 

L/L© 

HD19745*’ 

4700±100 

2.25±0.25 

1.5±0.2 

-15.4±1.0 

-0.05 

3.0+1.0 

470 

2.2+0.6 

1.9 

IRAS 12327-6523 

4200±100 

1.25±0.25 

2.0±0.2 

-18.9±1.0 

-0.40 

7.5+1.0 

310 

2+1.0 

2.9 

IRAS 13539-4153 

4300±100 

2.25±0.25 

2.0±0.2 

-87.3±1.3 

-0.13 

5.0+1.0 

630 

0.8+0.7 

1.6 

IRAS 17596-3952 

4600±100 

2.50±0.50 

1.9±0.2 

-27.2±2.3 

+0.10 

35.0± 5 

300 

1.7+0.6 

1.7 


^[Fe/H] = log e(Fe) - log e(Fe)Q with the solar iron abundance of log e(Fe) = 7.47 taken from Lodders (2003) 
^Castilho et al. (2000) obtained Tgff = 4750 K, log g= 2.5, = 1-2 and [Fe/H] = 0.0. 
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Table 2. Derived Li abundances and carbon istopic ratios. 


Abundance Feature 

IRAS 12327-6523 

IRAS 13539-4153 

IRAS 17596-3952 

HD19745 

LTE log e(Li) (Lil A 6103) 


4.15±0.3 


3.90±0.3 

NLTE log e(Li) (Lil A 6103) 


4.2 


4.0 

LTE log e(Li) (Lil A 6707) 

1.40±0.2 

4.10±0.3 

2.2±0.2 

3.70±0.3 

NLTE log e(Li) (Lil A 6707) 

1.6 

3.9 

2.3 

3.4 

12c/13c 

6.0±1 

>20 


16±2.0 






Relative Flux 


5 


Arcturus 



Wavelength (A) 


Fig. 1.— Spectra of Li-rich K giants and Arcturus near the Ho; line. Note the asymmetric 
Ha prohle of IRAS 12327-6523 






Fig. 2.— Observed (filled circles) and synthetic spectra (lines) of a region aronnd the Li I 
line at 6707.6 A for Arcturus, HD 19745, and IRAS 13539-4153. Note that the very strong 
Lil line at 6707.SA in HD19745 and IRAS13539 and absence of the same line in Arcturus. 
The stars have a similar effective temperature. The synthetic spectra with log e(Li) < —0.6 
and 3.7 provide the best fit for Arcturus and HD 19745, respectively. See text for discussion 
of synthetic spectra for IRAS 13539-4153. 






Wavelength (A) 


Fig. 3.— Observed (filled circles) and synthetic spectra (lines) of a region around the Li I 
line at 6707 A for IRAS 17596-3952 and IRAS 12327-6523. Spectra at 6707 A were predicted 
for the three different Li abundances as noted in the figure. The solid line shows the best fit 
to the observations: log e(Li) = 2.0 for IRAS 17596-3952 and 1.4 for IRAS 12327-6523. 







Wavelength (A) 

Fig. 4.— Observed (filled circles) and synthetic spectra (lines) of a region around the Li I 
line at 6103.5 A for Arcturus, HD 19745, and IRAS 13539-4153. Spectra at 6103 A were 
predicted for the three different Li abundances as noted in the figure. The solid line shows 
the best fit to the observations: log e(Li) < —0.6, 3.7, and 4.1 for Arcturus, HD 19745, and 
IRAS 13539-4153, respectively. 






8002 8004 8006 8008 

Wavelength (A) 

Fig. 5.— Determination of the ratio for three Li-rich K giants. In all the cases solid 

line is the best predicted spectrum for line at 8004.8 A resulting = 6, 16 and 

>20 for HD 19745, IRAS 12327—6523, and IRAS 13539—4153, respectively. 
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Fig. 6.— Li-rich giants in H-R diagram. Li-rich giants at the early-AGB phase (log L/Lq = 
2.8) and the luminosity function bump (log L/Lq = 1.8) from Charbonnel & Balachandran 
(2000) are shown by solid circles. The four Li-rich giants including HD 19745 from this study 
and one each from Reddy et al. (2002a) and Drake et al. (2002) are shown by open circles. 
Bigger symbols represent stars with log e (Li) > 3.3 and the smaller symbols represent stars 
with log e (Li) = 1.4 to 3.0. Evolutionary tracks of [Fe/H] = 0 (Girardi et al. 2000) for 
various masses are shown. 




